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Abstract 

We introduce the concept of Type-I/II generating functionals defined on the space 

of boundary data of a Lagrangian field theory. On the Lagrangian side, wc define an 
analogue of Jacobi's solution to the Hamilton-Jacobi equation for field theories, and we 
show that by taking variational derivatives of this functional, we obtain an isotropic 
submanifold of the space of Cauchy data, described by the so-called multisymplcctic 
form formula. We also define a Hamiltonian analogue of Jacobi's solution, and we show 
that this functional is a Type-II generating functional. We finish the paper by defining 
a similar framework of generating functions for discrete field theories, and we show 
that for the linear wave equation, we recover the multisymplectic conservation law of 
Bridges. 

1 Introduction 

The aim of the paper is to establish a theory of "generating functionals" for classical La- 
grangian field theories, playing a similar role to the generating functions in mechanics. 
Generating functions play an important role in mechanics, as a means of characterizing 
canonical transformations that preserve the structure of Hamilton's equations. A partic- 
ularly important generating function is the one associated with the time-dependent sym- 
plectic map that trivializes the dynamics of a Hamiltonian system, which is described by 
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the Hamilton-Jacobi equation. In turn, Jacobi's solution of the Hamilton-Jacobi equation 
can be viewed as a Type-I generating function which is given by the value of the action 
integral for paths that satisfy the boundary conditions and extremize the action. We will 
introduce a boundary Lagrangian for field theories that plays the role of Jacobi's solution, 
and demonstrate that functional derivatives of the boundary Lagrangian lead to expres- 
sions on the boundary of a space-time region that are analogous to the implicit expressions 
for a symplectic map expressed in terms of classical generating functions. 

This has important implications for discrete mechanics, since Jacobi's solution to the 
Hamilton-Jacobi equation is the exact discrete Lagrangian |18] , and the order of approxi- 
mation of a variational integrator can be characterized by the extent to which a computable 
discrete Lagrangian approximates the exact discrete Lagrangian. As such, a precise char- 
acterization of the boundary Lagrangian is essential for the construction and analysis of 
variational integrators for Lagrangian field theories. 

Geometry of the Space of Boundary Data To derive the expression for the bound- 
ary Lagrangian, we rely on the framework of Lawruk, Sniatycki and Tulczyjew [T3j, who 
describe the space of boundary data associated to a given PDE, showing among other 
things that a PDE determines an isotropic or sometimes even a Lagrangian submanifold 
of the space of all boundary data. Kijowski and Tulczyjew |12j connect these results with 
the polysymplectic description of field theories and give a description in terms of gen- 
erating functionals similar to the ones defined in this paper. In this formalism, a fixed 
Riemannian connection is needed to introduce the Hamiltonian formulation, but the intro- 
duction of metrics, connections, or other background objects might not be desirable when 
the goal is to describe a generally covariant theory (i.e., a theory which is invariant under 
arbitrary diffeomorphisms of spacetime; see |17j for a discussion of general covariance in 
field theories). One of the aims of this paper, therefore, is to recast these results in the 
multisymplectic formalism, where no extraneous background objects are needed for the 
Hamiltonian description. 

In recent years, Rovelli [22] has rediscovered many of these results in his search for a Hamil- 
tonian description of field theories in which no space-time split is made. He applied this 
formalism to general relativity using Ashtekar variables, rederiving among other things the 
Einstein-Hamilton-Jacobi equation. Our approach agrees with Rovelli's results wherever 
appropriate, but one of the advantages of the multisymplectic approach used here is that 
we are easily able to tie in our results with infinite-dimensional symplectic geometry. 

Inspired by these papers we introduce in this paper Type-I and Type-II generating function- 
als. While the technical details are given below, the idea is that for any given Lagrangian 
PDE we can define a functional Lgjj, and by taking the exterior derivative, a one- form 
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which takes Dirichlet boundary on dU (the elements of Vac/) into the corresponding Neu- 
mann data along dU (the elements of T*'^qij). We call the functional Lqjj appearing in 



(1.1) a Type-I generating functional. The image of this map turns out to be an isotropic 
submanifold of T*'^qjj (and in some cases even a Lagrangian submanifold) , so that the 
pull-back along dLgu of the canonical symplectic form uj on T*'^qij vanishes: 

{dLaufio = (1.2) 

If instead we prescribe mixed boundary conditions along dU , we refer to the correspond- 
ing functional as a Type-II generating functional. Taking variational derivatives of this 
functional will then yield the Dirichlet boundary conditions along the part of dU where 
Neumann boundary conditions have been specified, and vice versa. The resulting map 



again determines an isotropic submanifold of T*'^qu, so that a result similar to (1.2) holds. 



Generating Functionals and Multisymplecticity. The second goal of this paper is to 
shed further light on the multisymplectic form formula. This formula was introduced in [19] 
as a criterion for when the space of solutions of a PDE conserves a given multisymplectic 
form. Concretely, in the case of field theories derived from a Lagrangian function L, 
the multisymplectic form formula takes on the following form. If fi^ is the Lagrangian 



multisymplectic form (defined below in (1.5) and (1.6)), then for any subset U of the space 
of independent variables, and for any solution (p of the Euler-Lagrange equations defined 
on U , we have that 

/ (jV)* {ij^wh'v^L) = 0, (1.3) 

JdU 

where V and W are arbitrary first variations of (p. If we restrict to vertical first variations, 
this formula takes on the following form in coordinates: 

where V"", are the components of V and W , and the subscript v denotes differentiation 
with respect to x'^ . 

For Hamiltonian field theories on the other hand. Bridges f3] defined a different notion of 
multisymplecticity. In this context, a set of Hamiltonian PDEs is said to be multisym- 
plectic if they satisfy a certain differential conservation law. Under some restrictions, this 



conservation law can be rewritten in integral form to yield a result similar to (1.3), but the 
precise link between both formulations in full generality is not yet clear. 



In this paper we take a different approach to the derivation of (1.3) and its implications. 



Given a Lagrangian field theory, we construct the associated generating functional (1.1) 



and we relate this functional to the integral of the Poincare-Cartan form 0/, over the 
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boundary dU. In this way, we then show that the left-hand side of the isotropy condition 



(1.2), when written out in terms of the multisymplectic form 0.^ = —dQ^, is nothing but 



the multisymplectic form formula (1.3) 



The advantage of this approach is that it provides a criterion of multisymplecticity which 
can be applied to any set of PDEs for which a Dirichlet-to-Neumann map can be defined. 



From our point of view, therefore, the isotropy condition (1.2) is fundamental, and the 



multisymplectic form formula (1.3) arises as a consequence. 



Outlook: Discrete Lagrangian Field Theories. Our long-term goal is to clarify the 
concept of multisymplecticity for discrete Lagrangian field theories, in which the space of 
independent variables is replaced by a discrete mesh. While both the criterion of Bridges 
(see [U [5] and the references therein) as well as the multisymplectic form formulsj^ can be 
defined in this context, the relation between both is not clear yet. In the discrete context. 



we show that a generating functional akin to the one described in (1.1) can be introduced 
and that its image determines an isotropic submanifold of the space of (finite-dimensional) 
discrete boundary data. As in the case of continuous field theories, the condition of isotropy 
then gives rise to the discrete multisymplectic form formula. While we have not yet been 
able to establish the link between this condition and Bridges' discrete version of multisym- 
plecticity in full generality, we finish the paper with a simple discretization of the wave 
equation where this link can be established by direct computation. 



Dedication. We dedicate this paper to the memory of Jerrold E. Marsden. The methods 
pioneered by Jerry and his collaborators exerted a very profound influence on this paper: 
for the treatment of the infinite-dimensional geometry of the space of boundary data, 
we relied heavily on the foundational results from m [7], while the connection with field 
theory, and in particular the variational/multisymplectic formulation, uses the results from 
the GIMMSY manuscripts |10l 111] . In more recent years, Jerry was influential in the 
development of variational principles for discrete field theories, his paper [19j being the 
first to give a definition of multisymplecticity for a discrete Lagrangian field theory. 

The Geometry of Lagrangian Field Theories 

In this section, we briefiy recall the fiber bundle approach to classical field theory. We give 
a description in local coordinates; for a intrinsic description, as well as applications and a 
more in-depth discussion, see [H [6l Uni [H] and the references therein. 

Throughout this paper, we will consider fields as sections of a fiber bundle p : Y ^ X. 
Often, X will be spacetime and Y will be the product of X with a vector space V, but this 

^It is interesting to note from a historical point of view that a precursor of the multisymplectic form 
formula aheady appears in the seminal work of Courant, Friedrichs and Lewy [S]. 
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will not always be the case. We write the dimension of X as n+1, with n > 0, and we denote 
coordinates on X by {x^), = 0, . . . , n. We use the shorthand d'^^^x = dx^ A ■ ■ ■ A dx", 
and we put 

d^Xfj, := ig/Q^^.d^'^'^x 

= i-l^dx^ A • • • A dx^"-^ A A • • • A dx", (1.4) 

so that (up to sign) d^x^ is d^'^^x with the dx'^ term removed. 

On Y we choose coordinates (x'^, y"), a = 1, . . . ,k, adapted to the projection p, so that a 
section (j) : X ^ Y can be locally written as (p{x) = {x^, (j)"'{x)), with (/>"(x) locally defined 
component functions. 

We define the first jet bundle J^Y to consist of equivalence classes of local sections of 
vr, where two sections (p,(f)' : X ^ Y are equivalent at a point x S X if their first-order 
Taylor expansions around x agree, i.e., Tx(j) = Txtp' ■ On J^Y, we have local coordinates 
(x^ , y°' , y'^) , where the can be considered as the derivatives of with respect to the 
variables x'^, which we refer to as multi- velocities. Given a local section (p : X ^ Y p, 
we define the prolongation of (j) to be the section j^cj) of J^Y given in local coordinates by 

j'm = (xf^,nx),^ix) 

where the 0"(x) are the component functions of (p. 

By a Lagrangian density on J^Y we mean a map £ : Jiy ^ A"^^(^)> where A"^^(^) 
is the space of volume forms on X . Such a Lagrangian density can be written in local 
coordinates as £(x'^,y",y^) = L(x^, y", y^) (i""'"^x, where -L(x'^, y", y^) is a function on 
J^Y referred to as the Lagrangian function. 

Given a Lagrangian density £, we can introduce a number of geometric objects on the first 
jet bundle. In local coordinates, the Poincare-Cartan form is given by 

e, = (l - d^^'x + ^dy^ A d-x,. (1.5) 

We define the multisymplectic form Vt]^ on J^Y by 

= -dGi. (1.6) 



The expression (1.5) for the Poincare-Cartan form can be rewritten as 



31 

&L = Ld^-^'x + —{dy'' - y^dxn A d"x^, (1.7) 



since from (1.4) it follows that dx'^ A d^Xf^ = S'j^d'^'^^x, with 5^ the Kronecker delta. The 



advantage of this expression is that the forms dy"" — y^dx^ which appear in the second 
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term are contact forms, i.e., they vanish when pulled back along a prolongated section: if 
j^4>{x) = {x,(f)"'{x),d<p"'{x)/dx'^) in local coordinates, then 

- y^dx-) = d4>''{x) - ^^dx'' = 0. 
Along prolongated sections, we therefore have that 

(jV)*ei = (jVr(^d"+M- (1-8) 

This equality will often be useful later on. 



2 The Space of Boundary Data 

Let U be an open subset of X with boundary dU. We want to prescribe boundary data 
along dU with values in Y, the total space of the configuration bundle p : Y ^ X. We 
first describe the geometry of the space of all boundary data, and then we discuss some 
related spaces. We emphasize that at this stage, dU does not have to be a Cauchy surface, 
or even be spacelike: indeed, all of the definitions below are independent of the choice of a 
metric on X, so that in particular they can be applied to hyperbolic and elliptic problems 
alike. 

By an element of boundary data on U, we mean a section ip : dU -^Y of p, defined on 

the boundary of U. Wc denote by ^qu the space of all boundary data and wc now describe 
the tangent and cotangent bundles of this space. We can describe the tangent vectors S(p 
at a point G Vgt/ as follows: let be a curve in ^qu such that ife=o = and put 

dipe{x) 
d(p{x) := — 

for all X € dU. Note that the right-hand side takes vahies in the space of vertical tangent 
vectors: as each map ip^ is a section of p, we have that p o ip^ = Id, and therefore 

p{(Peix)) = 0. 

As the right-hand side takes values in we have that 6(p is a map from dU into VY 

with the property that dip{x) £ V^(^^^Y for all x G dU. These maps can alternatively be 
described as vector fields along or as sections of the pullback bundle (p*{VY). In any 
case, we have that the tangent space at a point is given by 

T^p'Sdu = {Sf -.dU ^VY : Sifix) G V^^^)Y for ah x e dU} . 



Tp{6<f{x)) = - 



6 



In local coordinates, an element Sip of T^'^qu is given by 

d 

where the 5ip"'{x) are locally defined component functions. Roughly speaking, we can think 
of the elements 5ip S T^^gjj as infinitesimal variations of the boundary data given by ip. 
Following [2], we will refer to T'^qu as the space of Cauchy data. 

We now describe the dual spaces T^'^qu. We restrict ourselves to the smooth dual, in 
other words, the space of smooth linear functionals from T^p^gu to M. Since the elements 
of T^^du are vertical vector fields along dU, the elements of the dual can be identified 
with linear maps vr from the space of vertical vector fields to the space of densities on 
the boundary, so that the duality pairing can be defined in terms of integration over the 
boundary by 

(599, vr) = / 7r-5(p. (2.1) 
JdU 

Here, tt ■ 6ip is the volume form on dU obtained by letting the linear map ir act on the 
vertical vector field 6ip. 

We now make this picture more precise. To describe T*ygi/, we consider first the tensor 
product bundle ip*{V*Y) (g) /\"^{dU), where V*Y is the dual bundle to the vertical bundle 
VY, Lp*{V*Y) is the pullback bundle over dU, and /\"{dU) is the bundle of volume forms 
(densities) on the boundary dU, viewed as a bundle over Y. Note that if there is a locally 
defined volume form dS on dU, then the elements of ip*{V*Y) (8) /\^{dU) can be written 
locally as TT = iTady^ dS. 

We now let T*'^gu be the space of sections of V*Y ® f\"'{dU), defined on the boundary 
dU. Each section vr G T*'^gij is a map from dU into ip*(y*Y) /\^{dU), with the property 
that 

^(x)Gi;%)y®A"(9c/) 

for all X £ dU. Again under the assumption that dS is a local volume form on dU, we 
have that tt can be written locally as 

tt{x) = TTady" dS. (2.2) 

For reasons that will become clear later on, we will refer to tt as the momentum in the 
direction normal to dU, and we will refer to the cotangent bundle T*'^gij as the space of 
normal momenta. 



The duality pairing (2.1) between T'^gu and T*^gu can then be written as 



{6(P,Tt) = / TT ■ S(p = 'Ka{x)5ip"'{x) dS. 
JdU JdU 
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For future reference, we point out that T*'^qi/ is equipped with a canonical weak symplectic 
form, given hy uj = —dQ, where is the one-form defined intrinsicahy as 

Q{7r){6Tr) = {TT{67r),TT) , (2.3) 

for all vr G T*'^qu and Jvr G Tt^{T*'^qu), where r : T*'^qu — ^qjj is the cotangent bundle 
projection (see [1]). Instead of this definition, we will often use the following defining 
property of 0: for every one-form (3 on ^qui we have that 

re = p, (2.4) 

where on the left-hand side /3 is interpreted as a map from into T*Vg^. 



3 The Boundary Lagrangian 

Let £ : J^Y A"^^(^) be a first-order La grangian density and denote, as before, 
^{i^4') = L[j^(f)) dV, with L(j^(p) a scalar function on J^l". To make the distinction 
between boundary data on dU and fields defined on the interior of U, we will denote the 
former by (/? G ^du^ while the latter are denoted by (p. 

We now define the boundary Lagrangian Lqjj as the functional on ^qjj given by 

Laui^) = I L{j'^)dV, (3.1) 



u 

where (/> is the unique section satisfying the Euler-Lagrange equations, and such that 
agrees with (p on dU : 

d ( dL , .1 A dL 
dxt' 

In other words, 



Lau{^) = S{4>), (3.3) 



where S is the action functional, and 99 and (j) are related as described in (3.2). This leads 
us to an alternative description of the boundary Lagrangian, which will often be useful in 
computations: as the solutions </> of the Euler-Lagrange equations are precisely the extremal 
points of the action functional 5, we see that Lgij{ip) is precisely the extremal value of S: 

Lqu{^)= ext 5(0), (3.4) 

<t>\dU=9 

where the extremal value of S is computed over the class of all sections (j) such that 
'P\dU — V'- Indeed, the variational characterization of the boundary Lagrangian is more 
broadly applicable than the characterization in terms of the action of the section that 
satisfies the boundary conditions and the Euler-Lagrange equations since in the case of 
degenerate Lagrangians, the section satisfying the boundary- value problem may not be 
unique, but the variational characterization remains well-defined. 
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The Space of Admissible Boundary Data. It will often happen that one cannot 
prescribe arbitrary boundary data for the Euler-Lagrange equations (3.2). We will denote 
by Xgjj the space of boundary data that give rise to a unique solution of the Euler-Lagrange 
equations, and we refer to %gu C ^qu as the space of admissible boundary data. The 
specification of %qu will depend on the type of PDE under study and the geometry of U. 
We give some examples at the end of this section. 



Functional Derivatives. We define the functional derivative of Lgu as follows: SLgu / Sip 
is the unique element of T*'^qjj such that 



dLguiip) ■ dip = [ ^^^-Sip 
JdU 



for every variation 5ip G T'^gu, where dLgu is the exterior derivative of the boundary 
Lagrangian. If not all boundary data are admissible, so that %gu is a proper subset of 
Vgu, then the admissible variations 6ip are elements of T%gu, so that dLgu{ip) G T*%gij. 

By applying the exterior differential to both sides of (3.3), we now obtain 

fi-Lgu{ip) ■ 6ip = dS{(p) ■ 6(1), 

where (p is the solution of the Euler-Lagrange equations with boundary data ip, and 6(f) is 
a first variation of (j), defined as follows. Let ipe be a curve in the space ^gu of boundary 
data such that ipe=o = f and -^fel^^Q = 6ip, and let (pe be the corresponding family of 
solutions of the Euler-Lagrange equations such that {(j)e)\gu = Ve- Then 6(1) is given by 

S<P[x) = — — , (3.5) 

Cl€ e=0 

for all X S dU. Explicitly, the first variations 6(f> satisfy the first-variation equation, 
obtained by linearizing the Euler-Lagrange equations. 

The exterior differential of the action functional is given by the first variation formula: 



dL ^ „ dL 



jn+1 



X 



dL d f dL\\ ^ f dL 



, , , <5y°d"+^2;+ / —6y^d''x.. 

u ydy" dx^'\^y^JJ Jgu dyf, 

As i;^> is a solution of the Euler-Lagrange equations, the integral over U vanishes and we 
conclude that 

f f)T 

dLgui^) ■6ip= —6y'' d^x^, (3.6) 
JdU dyl 
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and therefore 

^ = ^a^^y" ® ^* (d-x,) , (3.7) 

where l : dU X is the embedding of the boundary dU into X. Since (Px^ is an n-form 
on X , its puUback along i is a form of maximal degree on dU. If we choose coordinates on 
X that are adapted to dU, in the sense that dU is locally given by x^ = 0, then i is locally 
given by i{x^ , . . . , x") = {0, x^ , . . . , x"^), so that 

L* {(Txo) = (Txo, and i* = (i = l,...,n). 

As a result, in adapted coordinates we have that 

6Lgij dL 
dip dy^ 

A more intrinsic expression may be given when X is equipped with a metric tensor, as we 
now show. 



Normal Momenta. When a Riemannian or Lorentzian metric G on X is given, we may 
describe the functional derivatives as follows. In both cases, we have that 

L*{(rxf,) = dS, 

(see [23]), where is the outward normal to dU, dS is the induced metric volume form 



on dU, and indices are raised/lowered using the metric. The functional derivatives (3.7) 
can then be written as 

SLgu dL , „ 
TT := — - — = - — dy (g) do. 



6 if dy' 



so that by comparing with (2.2), we have for the components 

dL_ 



= ^'^M = Pant.. 



In other words, the boundary momentum Ha is the normal component of the spacetime 
momentum pa, so that we will refer to vr G T*'^qu as the normal momentum to the 
boundary dU . 



Multisymplectic Form Formula. A more careful derivation of (3.6) proceeds as fol- 
lows. The boundary Lagragian Lqu can be written in terms of the Poincare-Cartan form 
Gl as 



lau 
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where </) is the solution of the Euler-Lagrange equations with boundary data ip, and j^cj) is 



its first jet prolongation. Here, we have used (1.8) to bring in the Poincare-Cartan form. 



By taking the exterior derivative, we then obtain 

dlQuiv) -6^= [ {j'^r {ipyQi) , (3.8) 

JdU 



where the vector field F is a first variation of the solution (j), defined as before in (3.5). 
The advantage of this expression is that we can now take the exterior derivative again, to 
obtain 

(fLQui^) ■ {5^, V) = / (j V)* (ijiH/iiiy^^i) , (3-9) 

JdU 

where V and W are the first variations induced by 6(p and 5ip' , respectively. The proof of 
this results proceeds along similar lines as the proof of Lemma 5.1 in |llj . Since 
we now have that 

/ (jV)* (i,wi,ivf^L) = 0, (3.10) 

JdU 

for all solutions 4> of the Euler-Lagrange equations, and all first variations V, W. This is the 
multisymplectic form formula, first proposed in p!9]. Our derivation is close in spirit 
to the one in that paper, because of the link between the boundary Lagrangian and the 
action functional. In Section |4] we will see the multisymplectic form formula appear under 
a different guise, as the condition for the manifold of physical solutions to be an isotropic 
submanifold of the space of normal momenta. 

We finish by noting that our version of the multisymplectic form formula is somewhat less 
general than the one derived in ^19j, since we consider only vertical variations. However, 
as one of the motivations for this work is the derivation of a multisymplectic form formula 
for discrete field theories, for which infinitesimal horizontal variations are not really well- 
defined, this is not a fundamental restriction. 



Relation with the Crnkovic-Witten Symplectic Form. As pointed out by Rovelli 
[22], the Crnkovic-Witten symplectic form on the solution space can be related to the 
various structures on the space of boundary data. We present here a slightly different 
approach from Rovelli, emphasizing the link with the multisymplectic form formula. 

We assume that the base space X is equipped with a Lorentzian metric and we let U be 
a region bounded by two spacelike hypersurfaces, S+ and S_. We choose the orientation 
so that dU = — and we consider a boundary Lagrangian Lgjj for this particular 
geometry. In Wheeler's terminology [20J, U would be called a "thick sandwich." 



From the expression (4.1), to be proved below, for the pullback {dLgu)*uj, or alternatively 



from the multisymplectic form formula ( |3.10[ ) directly, we have that 

(iV)* (iiiH/ijiy^i) = / (j'V)* {ijiw'^jiv^L} 
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for all solutions (j) of the Euler-Lagrange equations, and first variations V, W . These ex- 
pressions, however, are nothing but the symplectic structures of Crnkovic-Witten [9] and 
Zuckerman [23] on the space of solutions associated to L, integrated respectively over 
E_ and E+. We have therefore shown that the definition of this symplectic structure is 
independent of the spatial hypersurface along which to integrate. This conclusion is of 
course not new, but the link with the multisymplectic form formula has hitherto not been 
established. 



Examples 



For most field theories, the boundary Lagrangian (3.1 ) cannot be computed explicitly. Here 
we present a number of examples where this is possible after all. 



Mechanics. In the case of a mechanical system on a configuration space Q and with 
regular Lagrangian L(g, v) : TQ — )• M, we let X be M and we take for Y the product M x Q. 
The projection p -.Y ^ X is then the projection onto the first factor. 

We let U be an open subset of X = M and we assume without loss of generality that U 
is an interval (0, h), so that dU = {0, h}. In this case, the space of boundary data is just 
the product Q x Q, where we think of the first, resp. the second factor, as specifying the 
configuration of the system at t = 0, h, respectively. Accordingly, the space of normal 
momenta is given by T*Q x T*Q. If L is nondegenerate and h is small enough, then it is 
well known that for every pair {qq, qi) in Q x Q there exists a unique solution q{t) of the 
Euler-Lagrange equations so that g(0) = qo and q{h) = qi (see, for instance, p^). 



According to the definition (3.1), the boundary Lagrangian is then given by 



Ldu{Qo,Qi) ■= [ L{q{t),Q{t))dt, 
Jo 

where q{t) is the unique solution of the Euler-Lagrange equations with boundary data 
(QOiQi)- 111 this way, we recover the exact discrete Lagrangian introduced in [18]. The 



variational derivatives (3.7) are 



and coincide with the momenta of the system at the begin and end point of the solution 
trajectory q{t), t G [0, h]. The minus sign is due to the orientation of the boundary {0, h}. 
Notice that these equations are precisely the implicit characterization of a symplectic map, 
where Lqu is viewed as a Type-I generating function. 
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Harmonic Functions. Consider secondly the case of harmonic functions on M", for 
which the Lagrangian is given by 

M0,V0) = ^||V0f 

so that the field equations are given by Laplace's equation, = 0. The action is then 
given by 

1 f 2 1 /■ d(t) 1 



^ Ju 2 Jqu dn 2 Ju 

where we have used the divergence theorem, with dcp/dn := n ■ V(/> the normal derivative. 
If we let (/> be a solution of Laplace's equation with prescribed boundary data (p on dU, we 
obtain for the boundary Lagrangian (3.1) 

Ldu{^) = l [ f^dl, (3.11) 
2 Jqu dn 



so that the variational derivative becomes 



d(j) (8) dl, 



6(p dn 

where dl is the Euclidian line element along dU. We conclude that the map which associates 
to each element of boundary data if the corresponding variational derivative 5Lqu/5(p is 
nothing but the Dirichlet-to-Neumann map of the Laplace equation. Last, if n > 3, the 
boundary Lagrangian (3.1) can also be formally related to the harmonic capacity of the 
domain U. 



The Wave Equation. Our last example concerns the wave equation (pu — 4>xx = 0, with 
Lagrangian 

m^,^.) = \{<f't-^l). (3.12) 

Our exposition follows |12| . We assume that {7 is a square in of unit length whose 
corner vertices are given by (0, 0), (1, 0), (1, 1), (0, 1) and we prescribe boundary data (^(t, x) 
along dU . An arbitrary solution to the linear wave equation can be written as (/)(t, x) = 
F{x — t) + G{x + t), where F and G are determined through the boundary conditions by 

F{x) + G{x) = ifiO, x), F{x - 1) + G{x + 1) = (/j(1, x), 

and 

F{-t) + G{t) = ip{t, 0), F(l -t) + G(l + t) = <f{t, 1). 

Note that for F and G to be determined from this set of equations, the boundary data 
needs to satisfy the following compatibility condition: 

(p{t, 0) + ip{l -t,l)- ip{0, x) - ip{l, 1 - x) = 0. (3.13) 
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For the unit square U, the space of admissible boundary data %du is therefore the space 
of ah smooth functions on the boundary dU that satisfy this condition. By substituting 
the solution of the wave equation with given boundary data if back into the action density 



(3.12), we then obtain the following expression for the boundary Lagrangian: 
Ldu{^) = / {^x{a, 0) - (/3t(0, a)) (99(1 - a, 1) - 99(0, a)) da. 





If the boundary of U consists (partly) of characteristic curves of the wave equation, then 
the wave equation will still be well-posed, but more restrictive compatibility conditions 
will arise. Assume, for instance, that U is the diamond shape in the (t, 2;)-plane bounded 
by the lines x it t = and x it t = 1. It will then only be possible to prescribe arbitrary 
boundary conditions along two adjacent sides of the diamond. In |12j . it is shown that 
the image of d.LQu{Lp) will in this case only be isotropic, while for the non-characteristic 
square the image of dLQu{ip) can be shown to be Lagrangian. We will not dwell on the 
continuous case any further and refer instead to [12^ , but later on we will see that a 
similar dichotomy arises in the discrete context. 



Maxwell's Equations. We finish with a brief note on the case of field theories with 
gauge freedom, and we treat electromagnetism as an example. The Lagrangian density for 
Maxwell's equations is given by 



C 



1 



(3.14) 



where F^^ is given in terms of the electromagnetic potential by F^y = d^Af^ — d^Ay, 
or in other words, F = dA when viewed as differential forms. As the Lagrangian does 
not involve the time derivatives of Aq, it is degenerate, and as a consequence Maxwell's 
equations, given by d^F>^'^ = 0, do not determine A^ uniquely. 

While this might preclude defining the boundary Lagrangian in a straightforward way as in 



(3.1), we can use the alternative definition (3.4) to introduce Lqu. Note first of all that any 
two solutions of Maxwell's equations A, A' differ by a gauge transformation. A' = A + df 



for some function /, while the action defined by the Lagrangian (3.14) is gauge invariant: 
S{A) = S{A -\- df). As a result, S takes on the same value for any solution of Maxwell's 



equation, so that we can use (3.4) to obtain Lqu by extremizing S. In differential form 
notation, we have that 



Lou {A) = / dAA*dA 
Ju 

= / d{A A *d^) - / AAd*dA 
Ju Ju 

A\Qu A {*dA)\Qjj, 
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where ^ is a solution of Maxwell's equations in the interior of U, so that d * dA = 0. On 
the last line, A\qu is the restriction of A to the boundary dU, and similarly for {*dA)\Qij. 

One issue now remains: how should the boundary data for A along dU be specified? The 
boundary Lagrangian inherits the gauge freedom of the original Lagrangian density: for all 
functions /, we have that LQu{A + df) = Lqij[A). Secondly, if ^ and A' arc boundary data 
whose tangential parts agree on dU, i.e. such that A^gu = A'^qu, then Lqu{A) = Lqu{A'). 
This can be seen by considering the following integral: since A^qjj — A'^gjj = 0, we have 
trivially that 

/ {A\9u-ALu)^{*d{A + A%Qu = 0. 

JdU 

By expanding the left-hand side we then obtain 

LguiA) - LguiA') - [ A' A i*dA)^9u + / A^Qu A (*d^')|9c/ = 0- 

JdU ' JdU 

Through repeated use of Stokes' theorem and the Maxwell equations, the last two terms on 
the right hand side can be made to vanish, so that we can conclude that Lgu{A) = Lgu{A'). 

As a result, rather than specifying the four-dimensional vector potential A = A^j^dx^ along 
dU, it is sufficient to specify only the tangential part A^qij. Moreover, since Lqjj is gauge 
invariant, Lgu depends only on the equivalence class of A^qu under gauge transformations 
and so we may define a new functional 

Ldu{dA\Qij) := Ldu{A\gu), 

where we assume that the topology of dU is sufficiently simple, so that we can identify 
exact two-forms with equivalence classes of one-forms up to exact forms. If dU is locally 
given by x° = 0, then dA\g^ can be identified with the magnetic field B. In general though, 
dA^gjj will contain both electric and magnetic parts relative to the splitting of Minkowski 
space into space and time, but we conclude that it is sufficient to specify B := dA^gjj along 
dU. 



4 Type-I/II Generating Functionals 

The purpose of this section is to show that the boundary Lagrangian Lgjj can be viewed as 
a generating functional, in some appropriate sense. We proceed along the following route: 
first, we take variational derivatives to define a mapping from ^gjj to T*']^gu, given by 

dLgu :ipe%u^ ^ 

and generated by Lgjj. That is, dLgif{ip) gives us the normal momenta 6Lgu/5Lp, given the 
fields if at the boundary. Secondly, we show that the image of this map is a Lagrangian 
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submanifold of T*'^qij with respect to the canonical symplectic structure (2.3). Since 



dLgij is an exact one-form, this will be a straightforward consequence of the property 



(2.4). Based on these two properties, we will say that Lgu is an example of a Type-I 



generating functional. As an interesting side-result, we obtain that the Lagrangian nature 



of the image of dLgu is intimately related to the multisymplectic form formula (3.10). 

In the second part of this section, we show that the image of dLgu can be generated by 
other functionals as well. More precisely, we can imagine subdividing the boundary dU 
into two subsets A and B, and prescribing boundary data (pA along A and normal momenta 
ttb along B. We can then define a functional Hqu{(pa-, t^b) of these data, whose variational 
derivative with respect to ipA is the normal momentum along A, and whose variational 
derivative with respect to ttb is the field (p along B. Consequently, the image of d(3n 
coincides with the image of dHgu in T*'^qu, and in analogy with mechanics, we will call 
Hqu a Type-II generating functional. 



Geometry of Generating Functionals. Let & : ^qjj — )• M be an arbitrary functional 
on the space of boundary data. The exterior derivative d(3 is a closed form, and its image 
is a Lagrangian submanifold of T*ygu. If instead we restrict d& to a subspace %du of 
T*'^Qij, the resulting image is no longer Lagrangian, and we obtain merely an isotropic 
submanifold. We collect these observations in the following theorem, the proof of which is 
a combination of standard results in symplectic geometry (see, for instance, [T]). 

Theorem 4.1. Let %dU C ^su be a subset of the space of boundary data, let & be a 
functional on OCdjj, and consider the image "Mqu '■= d&{3CQu)- Then l^lgjj is an isotropic 
submanifold ofT*'^Qu. If JCgu = ^qjj, then JAqjj is Lagrangian. 

By analogy with mechanics, we refer to any functional S on the space of boundary data 
as a Type-I generating functional. We now make this analogy more precise. Consider 
a mechanical system with configuration space Q, and identify the space of boundary data 
with QxQ and the space of normal momenta with T*QxT*Q, equipped with the symplectic 
form — r^o © ^ii where Oj, i = 0,1, is the canonical symplectic form on the ith factor. A 
function S{qo,qi) on Q x Q then generates a symplectic transformation 

dS\ / dS 



dqoj V 9qi 

in the standard sense, and our definition is an extension of this concept to classical field 
theories. 



The Boundary Lagrangian as a Type-I Generating Functional. Clearly, the 
boundary Lagrangian Lqu : ^qu — )• M is a Type-I generating function in the sense de- 
scribed above. We let J^qjj be the image of dLgu, restricted to the space of admissible 
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boundary conditions OCdU- We can identify 'Mqu with %qu, and under this identification 
the restriction to 'Mou of the symplectic form u on T*'^qu is given by the puh-back form 



{dLQu)*uj. From Theorem 4.1 we know that J^dU is isotropic, so that this form vanishes. 
However, some interesting results can be obtained by exphcitly writing out (dLg^)*a; and 
equating the result with zero. 

The canonical symplectic form oj is given by a; = — dB, where is the canonical one- form 
on T*'^Qu defined in (2.3). The pullback of along dLgu : ^qu — )• T*'^qu then satisfies 
{dLQjj)*Q = dLgu because of (|2.4[), so that 



((dLac/)*e) (^) • 6ip = dLdui^) ■Sip= [ {j'^r (iiiyQi) 

JdU 



where we have used (3.8). Here, (f) is again the unique solution of the Euler-Lagrange 
equations with boundary data ip, and V is a first variation of (p induced by the boundary 
variation 6(p. By taking exterior derivatives of both sides and using (|3.9l), we obtain 



{{dL9uruj){ip)-{Sip,6ip') 



au 



(4.1) 



but by isotropy, we have that the left-hand side of this formula has to vanish. As a result. 



we conclude that 'Mgjj is isotropic if and only if the multisymplectic form formula (3.10) 
holds. 



The De Donder-Weyl Equations. Given a Lagrangian density C = L{x^,y'^,y'^) dV, 
we define the multi- momenta pa (where a = l,...,k, ^ = 0,...,n) and a scalar momentum 
p as 

dL dL 
= — and p = L- ^^y",. (4.2) 



These momenta can be defined intrinsically by considering the Legendre transformation as 
a map from the jet bundle to its extended dual (see, for instance, |6j). We now introduce 
the Hamiltonian function T-L^x^ jy"" ,p,pa) as 

^(x^y^p,P^^) = ext [p + p^aV"^ - L(x^y^y^)] , 

where we take the extremum over all values of y^. For a related, coordinate-invariant 
definition oi Ti, we refer to j^. Note that 7i vanishes identically on the image of the 



Legendre transformation (4.2), and that L does not necessarily have to hyperregular for Ti 
to be defined. 

For the remainder of this paper, we will focus on the locally defined function H{x'^, y"',Pa) = 
Pa.y'^ — i(x^,y",y^), which we will term the multi- Hamiltonian (function). We em- 
phasize, however, that unlike Ti, H is in general only locally defined, but all the objects 
defined in this section can be defined in terms of % solely. 
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We now rewrite the action density S in terms of the Hamiltonian function: 

Ju 

and we notice that by taking variations of y"" and pa, we obtain 

nsMy',^J . ^ X (- - i^) + - i) «) ''-^ 

+ / p^^Sy'^rx^. (4.3) 

Under the condition that the variation Sy"" vanish on the boundary dU, we obtain the 
fohowing set of partial differential equations (referred to as the De Donder- Weyl equa- 
tions) : 

dr^dH_ ^^^^ dpl^_cm 

dxi^ dpa ' dxi^ dy°- 

Type-II Generating Functionals. We consider again a fixed domain U C X and we 
divide the boundary dU into two disjoint parts A and B: dU = AUB. We suppose that, on 
A, we are given fixed boundary fields ipA, while on B we are given fixed normal momenta 
ttb- We recall that the components of the normal momenta can be expressed as 

(7^B)a=P^*(^^X)> 
where l : dU ^ X is the inclusion of the boundary in X. 

For given boundary data {(fA, t^b), let {(p'^{x),pa{x)) be the solution of the De Donder- Weyl 



equations (4.4) with those boundary data, and define the functional 

H9u{ipA,TTB) = -SH{(l)^,p'i) + / (vTb) • 

Jb 

= - [ {Pa^ct>':^-H{r,pj:))d^-^^x+ j p/rd^x^, 

Ju Jb 

which we refer to as the boundary Hamiltonian. We now compute the derivative of Hgu, 
keeping in mind the boundary conditions, so that 5ip^A = '^'^IB = 0. A similar computation 
as for the derivation of (4.3) yields 

DHM.A, • iS,A, S.B) ^ I (- + Sy^ + (|g - ^) V.') c^-^. 

+ / Spi^y'^d'^x^- [ jfjy'^d^x^. 
Jb J a 



18 



The integral over the interior vanishes since {(p"',pa) is a solution of the De Bonder- Weyl 
equations, and the boundary integrals can be written in terms of the normal momenta as 

DHgu{ipA,'KB) ■ {S^PA^^TTb) = / StTb ■ 4>\B - / T^\A ' ^^A, 

Jb J a 

so that the variational derivatives are given by 

— = -7r\A, and — = (p\B- (4.5) 

We compare this with the case of mechanical systems. Given a finite time interval [to,ii], 
a Type-II generating function is a function S{qo,pi) depending on the position variables 
qo at the initial time Iq and on the momenta pi at the final time ti . The final position qi 
and the initial momentum po are then defined by 

as , as 

Po = and qi = ——. 

oqo opi 



The relations (4.5) are the analogue of these expressions for field theory, where the relative 



minus sign is again due to the orientation of the boundary. 



5 Example: Euler Discretization of the Wave Equation 

In this section, we show that concepts such as the boundary Lagrangian can also be defined 
for discrete field theories. Under a few modest assumptions on the discretizations used, we 
recover, among other things, the discrete multisymplectic form formula of ^19j. Throughout 
this section, we use the Euler discretization of the linear wave Lagrangian as a motivating 
example. For this special case, we discuss the influence of the geometry of the boundary on 
the expression for the multisymplectic form formula, and we show that the latter reduces 
to the multisymplectic conservation law derived by Bridges et al. jSHS]- Our treatment 
is inspired by the one in [13\ I19j . where different discretizations and more complex field 
theories are treated. 

For the sake of convenience, we restrict ourselves to a scalar field theory with two inde- 
pendent variables, which we label by t and x. We assume that we are given a regular 
quadrangular mesh in the base space, with mesh lengths At and Ax, and we denote the 
nodes in this mesh by (n, i) G Z x Z. Note that the node (n, i) corresponds to the point 
:= (nAt, iAx) in E?. We denote the value of the field u at the node (n, i) by uf. 

As in |19| , we introduce a discrete version of the jet bundle as follows. We define a triangle 
at (n, i) to be the ordered triple ((n, i), (n,i + 1), (n + which we denote by Af (or 

simply by A if no confusion can arise), and we let X^^ be the set of all such triangles. 
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Given a triangle A = ((n, z), (n, i + 1), (n + 1, z)), we refer to the vertices in a concise way 
by Ai := (n, i), A2 := (n, i + 1), and A3 := (n + 1, i). 

The discrete jet bundle is then given by (see 



JiY := {«, € M3 : ((n, i), (n, i + 1), (n + G Xa}, 

and so is equal to x M^. Given an element (n", u2_^_i,u^^^) of the discrete jet bundle, we 
define the triangle midpoint by u" := {u^ + u^_^i + u^~^^)/3, and we introduce the following 
first-order expressions for the temporal velocity and the spatial velocity w": 

„,n+l _ „.n „.n _ n 

z;r:= - , ' and < := !^i±l-^. (5.1) 
Consequently, we can discretize a given Lagrangian density L on J^y by 

L,(u,,n,+i,tx, - , , ^ • (5.2) 



For the linear wave equation, with Lagrangian (3.12), the discrete Lagrangian becomes 



L,«,<,i,nr^) 



,^.1, _ AtAx //n^^V 



At J V Ax 



Poincare-Cartan Forms. Given a discrete Lagrangian L^, we introduce the discrete 
Poincare-Cartan forms 0]^, 0| and by 

ei«,nr+i,<+i) := DiL,(nr,<+i,<+i)d<, 

and similarly for G| and B^. Note that these forms are one-forms on the discrete jet 
bundle, and that el + el + Ql = dL. Furthermore, we put n'l = -dQ^ (for A; = 1, 2, 3), 
so that 

nl + nl + nl = o. (5.3) 

For the linear wave equation, a straightforward computation yields 

ni = ^((Ax) dv^ A d< - (At) dw^ A d<) (5.4) 



as well as 



nl = ^{At)dw^ A dul, and J^i = -^(Ax) dff A d<, (5.5) 



where vf and are given by (5.1). 
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Figure 1: Three adjacent triangles touching a common vertex, labeled by n". 



Discrete Euler-Lagrange Equations. Given a finite subset U C X/^ of the space of 
triangles, we form the discrete action sum as 

Su{u) := ^ L,(j^n(A)). 
Aec/ 

Here, n is a discrete fields assigning to every node (n, «) a field value u", and j^u is its 
first jet extension, defined by 

(jin)(A) = KAi),7x(A2),n(A3)) G 4^. 

We now focus on a particular configuration [/, consisting of three adjacent triangles, as in 
Figure [Tj The action sum for this U is explicitly given by 

(5.6) 



Su{u) = L,{u^, ^) + <, + L,«-\ <) 



,n— 1 „,n—l „,n\ 



By keeping the values of the field on the boundary fixed, and taking variations with respect 
to Uni, we obtain the following discrete Euler-Lagrange equations: 



0. 



, "-j+i , "-i ; — ^- (5.7) 
We can rewrite these equations in terms of the discrete Poincare-Cartan forms as 

ei(A?) + ei(Ati) + ei(Ari) = o, (5.8) 

where A", etc., refer to the triangles defined in Figure [Tj 



For the wave equation (3.12), the discrete Euler-Lagrange equations result in the standard 
second-order scheme 



n-l 



+ 



<+i-2< + <_i 
(Ax)2 



0. 



(5.9) 
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Discrete Boundary Lagrangian. We now mimic the construction in Section [3] to in- 
troduce a discrete version of the boundary Lagrangian. We refer again to Figure [T] and we 
define the boundary Lagrangian to be the extremal value of the action sum ( |5.6[ ), taking 
variations over the interior node m^. For the sake of notation, we denote the values of the 
field on the boundary by uqu, so that 

,,„^r — („.n n+l n+l n n-1 „.n-l\ 

The boundary Lagrangian is then given by 

^S^/(^aL^):=ext[L,(n^<+l,<+l)+L,(ntl,<,<_Y) + ^d(<"^<+l^<)], (5.10) 
where the extremum is taken over all uf. Alternatively, Lgu can be defined by solving the 



discrete Euler-Lagrange equations (5.7) for given the values of uqu as boundary data. 



and substituting this value into the action sum (5.6). 



Discrete Multisymplectic Form Formula. We now derive the discrete multisymplec- 
tic form formula by twice taking the exterior derivative of the boundary Lagrangian Lqjj. 
By taking the exterior derivative of both sides of (5.10) and using the definition of the 
Poincare-Cartan forms given above, we obtain 

k=l 1=1 



Using the discrete Euler-Lagrange equations (5.8), we can rewrite this as 

k=l l=l;l=/=k 

and by taking another exterior derivative of both sides, and using that d^L, 



dU 



finally obtain 







k=l , 



3 

E 

= l;lj^k 



(5.11) 



0, we 



(5.12) 



This is precisely the multisymplectic form formula derived in [19] , applied to the triangular 
domain of Figure [!} For the linear wave equation, we can substitute the expressions (5.4) 
and (5.5) for the discrete multisymplectic forms to obtain 

-^{dWi A d< - dwti A duti) - -^^idv^ A d< - dv^-^ A du^-^) = 0. 

This is the multisymplectic conservation law of Bridges and Reich ^ applied to the linear 
wave equation, but we caution against taking this result too far: in general the discrete 
Euler-Lagrange equations (5.7) will be different from the discrete multisymplectic equations 
of Bridges and Reich. Consequently, the discrete multisymplectic form formula (5.12) will 
also lead to different multisymplectic conservation laws. 
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Lagrangian Submanifolds and Characteristics. With the notations of the previous 



paragraph, let JAgjj be the image of dLgu given by (5.11 ). Note that 'Mqu is a submanifold 
of (T*Q)^®, where we have one copy of T*Q for each boundary point in Figure [l] Note 
that JvIqu is the image of an exact one-form, and hence determines at least an isotropic 
submanifold of {T*Q)^^. If there are no extra conditions on the boundary data, then JAqu 
will be Lagrangian. 

We now discuss the relation between the location of the boundary data and the nature of 
Mac/. To keep the computations to a minimum, we restrict ourselves to the linear wave 



equation. We first solve the discrete wave equations (5.9) to determine u'^ in terms of the 



boundary data. The discrete wave equation can be rewritten as 

2{c' - !)< = c2«+i + - + <-^) (5.13) 

where c := At/ Ax is the aspect ratio of the mesh. If the CFL condition is satisfied, so 
that c < 1, this expression can be used to determine tif in terms of the boundary data. 
However, when c = 1, the left-hand side vanishes and we merely obtain a relation between 
the boundary data. In this case, the boundary data are located on characteristics of the 
continuous wave equation. 

As a result, whenever the boundary data is characteristic, there exist supplementary com- 
patibility conditions between the boundary data, and hence 'Mqjj is strictly isotropic. When 
the boundary is noncharacteristic, there are no further conditions on the data and M is 
Lagrangian. We outlined a similar result for the continuous wave equation at the end of 
Section |3l 



6 Conclusion and Future Directions 

In this paper, we consider the space of boundary data for a Lagrangian field theory, whose 
tangent space is the space of Cauchy data, and we introduce a duality pairing between the 
space of Cauchy data and normal momenta. We then introduce the boundary Lagrangian, 
which is the analogue of Jacobi's solution of the Hamilton-Jacobi equation, and consider 
more generally the concept of Type-I and Type-II generating functionals for a Lagrangian 
field theory. 

We are interested in the following topics for future work: 

• Continuous and Discrete Hamilton-Jacobi Theory for Field Theories. By comput- 
ing variations of Jacobi's solution of the Hamilton-Jacobi equation, one can recover 
the Hamilton-Jacobi equation. In a similar way, we intend to use the concept of 
the boundary Lagrangian to systematically derive a Hamilton-Jacobi equation for 
covariant field theories. 
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In [2T], a discrete Jacobi's solution was used to derive a discrete Hamilton-Jacobi 
equation, which is in turn related to the Hamilton-Jacobi-Bellman equation of op- 
timal control. It would be interesting to develop the analogous connection between 
discrete Hamilton-Jacobi theory for fields and numerical methods for the optimal 
boundary control of Lagrangian field theories. 

• Variational Error Analysis of Variational Integrators for Field Theories. Since the 
boundary Lagrangian Lqu and the boundary Hamiltonian Hgjj are exact generating 
functionals for Lagrangian and Hamiltonian field theories, it would be natural to 
extend the theory of variational error analysis introduced in [18] to the setting of 
discrete field theories, and to develop general techniques for constructing variational 
integrators for field theories by extending the approaches for constructing Lagrangian 
and Hamiltonian variational integrators described in |15| [T6] . 

• Connections with Multisymplectic Integrators for Hamiltonian PDEs. A variational 
characterization of discrete multisymplectic variational integrators for Hamiltonian 
field theories is a natural direction to pursue, and will involve a combination of the 
insights developed in this paper, and a generalization of the methods developed in 
[T6] . The goal would be to provide a systematic characterization of the associated 
discrete multisymplectic conservation laws, and appropriate invariance properties of 
the discrete Hamiltonian which would lead to a discrete multisymplectic Noether's 
theorem. This connections between such a discrete variational approach and the 
multisymplectic integrators described in [H [5] remain to be elucidated. 
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